A junction of three quantum wires: restoring time-reversal symmetry by interaction 
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We investigate transport of correlated fermions through a junction of three one-dimensional quan- 
tum wires pierced by a magnetic flux. We determine the flow of the conductance as a function of 
a low-energy cutoff in the entire parameter space. For attractive interactions and generic flux the 
fixed point with maximal asymmetry of the conductance is the stable one, as conjectured recently. 
For repulsive interactions and arbitrary flux we find a line of stable fixed points with vanishing 
conductance as well as stable fixed points with symmetric conductance (4/9)(e 2 //i). 
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Electronic transport through quasi one-dimensional 
(Id) systems is of current experimental and theoretical 
interest. In one spatial dimension correlations play an 
important role and the physical properties of interact- 
ing fermions show distinctive non-Fermi liquid features. 
Generically such systems can be described as Tomonaga- 
Luttingcr liquids (TLLs) characterized by a vanishing 
quasi-particle weight and power-law scaling of correla- 
tion functions. pj For spinless fermions, on which we fo- 
cus here, the characteristic exponents are expressible in 
terms of the interaction dependent TLL parameter K. 
For repulsive interactions < K < 1, while K > 1 in 
the attractive case. In TLLs inhomogeneities have a dra- 
matic effect as can be inferred from the singular behavior 
of response functions of homogeneous models. 0,0 

In an important first step the conductance G through 
a TLL with a single impurity was understood. 0, 0. For 
vanishing energy scale s (e.g. temperature) and < K < 
1, G tends towards following a power-law. The low- 
energy physics is governed by the "decoupled chain" fixed 
point (FP) . The scaling dimension of a hopping term con- 
necting two open ends is 1/K and leads to G oc s 2 ^ 1 '-^ -1 ). 
One can understand this behavior in a simple picture. 
Due to the interaction the self-energy develops long range 
oscillatory behavior and the scattering off the resulting 
effective impuritypotential leads to the power-law sup- 
pression of G.[fj|, lD, la] For K > 1, the conductance ap- 
proaches the impurity free limit. Close to the "perfect 
chain" FP an impurity has scaling dimension K and 
the correction to the impurity free conductance scales 
as s 2 ^^ 1 ). In this case the effective impurity potential 
leads to a resonance at the chemical potential. 

Recently junctions of several quasi Id quantum wires 
were realized experimentally with single-walled carbon 
nanotubes. 0, ^} They might form the basis of electronic 
devices. Already the physics of the three wire junction 
(Y-j unction) is considerably richer than the one of a sin- 
gle impurity. Taking into account correlations transport 
through such systems was investigated theoretically in 
Refs.llll andll2l These studies posed a number of inter- 
esting questions. In Ref. ^3 a symmetric triangular Y- 
junction pierced by a magnetic flux </) (measured in units 
of the flux quantum hc/e) was studied. In this geometry 
time reversal symmetry is broken. In the non-interacting 



case this generically leads to an asymmetry of the conduc- 
tance from wire v to wire v' and vice versa: G v ^> ^ G v > 
For K > 1 at flux <f> — ±ir/2 one FP was found applying 
an exact method adopted from boundary conformal field 
theory. The FP corresponds to the case of maximal asym- 
metry of G u y. We here consider TLL wires connected to 
semi-infinite Fermi liquid (FL) leads with TLL-FL con- 
tacts that are modeled to be free of fermion backscat- 
tering. The results of Ref. ^3 obtained for semi-infinite 
TLLs can be extended to our modeling. Then maximal 
asymmetry is given by G u y = 1(0) and G,/ „ = 0(1) in 
units of the conductance quantum e 2 //i.[L| The scaling 
dimension of the most relevant operator at this FP is 
A = 4K/ (3 + K 2 ) and the correction to the FP conduc- 
tance scales as s 2 ^" 1 )^ For 1< K < 3, 2 (A - 1) > 
and the "maximal asymmetry" FP is attractive. This 
implies that independently of the junction parameters at 
low energy scales the conductance is 1 from v to v 1 and 

from v' to v or vice versa. It was conjectured that for 

1 < K < 3 and all (/>, except \<t>\/ir being integer, this FP 
is the only stable one. 

We investigate the same physical problem considering 
arbitrary fluxes and attractive as well as repulsive inter- 
actions. An approximation scheme based on the func- 
tional renormalization group (fRG) method is used. It 
was earlier applied successfully to transport problems in 
inhomogeneous TLLs. By comparison with numerical re- 
sults and exact scaling relations this scheme was shown 
to be reliable for 1/2 < K < 3/2. In particular, the 
scaling dimensions of the two FPs of the single impu- 
rity problem come out correctly to leading order in the 
interaction J7.0,@] For the Y-junction we here confirm 
the conjectured stability of the "maximal asymmetry" 
FP for all (j>, with \4>\/ir not integer, and reproduce the 
scaling dimension A to leading order in U. For U > 
and arbitrary <j> a line of stable FPs with G = is iden- 
tified. In cases with symmetric conductance we suppress 
the indices on G. We surprisingly find additional "per- 
fect junction" FPs with symmetric conductance G = 4/9 
that for U > are stable and have a scaling dimension 
not discussed before. In a Y-junction of non-interacting 
wires without flux G = 4/9 is the value maximally al- 
lowed by symmetry. Although time reversal symme- 
try is explicitly broken, for systems that flow towards the 
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"perfect junction" FPs the electron correlations cause the 
conductance - an observable that is commonly believed 
to indicate this breaking of symmetry - to behave as if 
the symmetry is restored. In a certain sense (see below) 
this also holds for systems that flow towards the line of 
G = FPs and thus for all parameters. 
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FIG. 1: Symmetric junction of three quantum wires. 



Our wire Hamiltonian is 



H v = (ct ^+i^ + h.c.J 
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in standard second-quantized notation. The different 
wires are indicated by an index v = 1,2,3. The am- 
plitude of the nearest neighbor hopping is set to 1. The 
nearest neighbor interaction Uj is allowed to depend on 
the position. It is set to zero for j > N, i.e. the wires of 
length N are connected to non-interacting leads. Close 
to the contacts the interaction is switched off (spatially) 
smoothly to avoid any fermion backscattering.J^ The 
bulk value of the interaction is denoted by U. We here 
consider the half-filled band case. To prevent depletion 
of the interacting region we shifted the operator rij v by 
the average density 1/2. The model with interaction U 
across all bonds (not only the ones within [1, AT]) shows 
TLL behavior for \U\ < 2 withQ 
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The Y-junction sketched in Fig. is described by 



H Y = -ty ^ (c\ v cq jV + h.c.j + V np lV 

3 

-t A ]T (>/ 3 c^c ,, + i + h.c.) , (3) 

where we identify the wire indices 4 and 1. 

Our starting point to calculate the linear response con- 
ductance of Ylv=i Hi> + Hy is an exact hierarchy of dif- 
ferential flow equations for the self-energy E A and higher 
order vertex functions (in the presence of the interac- 
tion, the leads, and the junction), with an infrared en- 
ergy cutoff A as the flow parameter. It is derived us- 
ing the fRG. The hierarchy is truncated by neglecting 



n-particle vertices with n > 2, and the 2-particle ver- 
tex is parametrized by a renormalized nearest neighbor 
interaction U A . This implies that terms of order U 2 
are only partly taken into account and that E A is fre- 
quency independent. The important spatial dependence 
of E A , is however fully kept. The details of this proce- 
dure are given in Refs.Q and^. At the end of the flow (at 
A = 0), the self-energy can be regarded as an effective, 
A-dependent impurity potential E with non-vanishing 
matrix elements Ejj and Ejj+i, where j is restricted to 
the interacting region. Due to the symmetry of the junc- 
tion the matrix elements are the same for the three TLL 
wires. We here focus on the zero temperature case for 
which the flow equations can numerically be solved for 
up to N = 10 7 lattice sites. 0, Generically Y>j j and 
Ej J+ i oscillate around an average value with an ampli- 
tude that decays slowly with increasing distance from the 
junction. 

The conductance G„,„' = \t^^'\ 2 can be calculated 
from the effective transmission t vv > (at the chemical po- 
tential /i = 0), which in turn can be expressed in terms 
of real space matrix elements of one-particle Green func- 
tions of the system. Using single particle scattering 
theory, |8( the conductance can be written as 



4(Im 5 r 



\g 3 — 3g + 2cos< 



(4) 



with g = (— V — ty Qi \)/\t^\. The Green function Q is 
obtained by considering E as an effective potential for a 
single wire setting ty = and Qi i denotes its diagonal 
matrix element taken at site j '• = 1 . It is evaluated at en- 
ergy e + iO with e — > 0. Eq. Q holds if v, v' are in cyclic 
order. G v i tV follows by replacing 4> — > —<f>. For symmetry 
reasons we only have to consider O<0<7r/2. In par- 
ticular Eq. Q can be applied for U — and shows that 
the conductance can be parameterized by the flux and 
a single complex parameter g. Via the flow of the self- 
energy Gi,i for U ^ develops an additional dependence 
on (tyjt^, V), U, and most importantly on N. The en- 
ergy scale S n — irvp/N (with the Fermi velocity vp) is 
a natural infrared cutoff of our problem. || To obtain a 
comprehensive picture of the low-energy physics we in- 
vestigate the flow of g as a function of Sn and use Eq. 
Pfl to calculate G v y for a given g. In Fig. the flow of 
g is shown for different <f> (and U = —1). Each line is ob- 
tained for a fixed set of junction parameters (ty, t&, V) 
with N as a variable. As Img has the opposite sign of 
Im(/i,i < it is restricted to positive values. 

Eq. (0J allows for three special conductance situations 
which turn out to be the FPs of the flow. On the real 
axis (green in Fig. [5J , the conductance G vanishes except 
at special, ^-dependent points. They are given by the 
crossings of the real axis with the set C{<j>) (red in Fig.EJ, 
on which the reflection R = 1 — \t v y\ 2 — \t u >. v \ 2 takes a 
local minimum. For ip — 0, C(0) is given by the circle 
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(Reg + 1/2) 2 + (Im 5 ) 2 = (3/2) 2 . For > 0, C(<f>) has a 
more complex analytical form not presented here. At the 
crossings the conductance is symmetric with G = 4/9. 
They are located at g = —2 cos (0/3), 2 cos ([it ± </>]/3) 
and in Fig. [2] are indicated as blue circles and a yellow 
diamond. As a peculiarity of <f> = on C(0) one finds 
G = 4/9 for all Img, not only Img = 0. A situation with 
G = is also reached for \g\ — > oo. For < < 1/2 
and 5 = e** Eq. gj) yields = 1 and G„> jV = 0. At 

this point, indicated by cyan squares in Fig. [2 maximal 
asymmetry of the conductance is achieved. 

We next discuss the FP scenarios depicted in Figs. 0a)- 
d). For U ^ the general form of the flow diagrams 
is independent of the absolute value and sign of U. In 
Figs. [21 a)-d) results for U = — 1 are shown. At <f> — 
we find two "perfect junction" FPs with G = 4/9 and a 
line of "decoupled chain" FPs with G = (green line). 
For all U < the "perfect junction" FP at g = 1 (yel- 
low diamond) is the only stable FP. All trajectories are 
attracted towards C(0) and reach this FP following C(0). 
For all U > it turns unstable and the line of "decoupled 
chain" FPs is stable. In addition the "perfect junction" 
FP at g = —2 (blue circle) is stable. The basin of attrac- 
tion of g = —2 is given by C(0). Increasing (f> from the 
"perfect junction" FP at g — 1 splits up into three FPs 
- the two "perfect junction" FPs at g = 2 cos ([ir ± </>]/3) 
and the "maximal asymmetry" FP at g = e % ^ . For all 



U < the latter is the or% stable FP but becomes unsta- 
ble for U > 0. A third "perfect junction" FP is located 
at g = —2 cos (0/3). For all U > each of the three 
"perfect junction" FPs has a basin of attraction given by 
one of the three parts of C(cf>) which are separated by the 
"maximal asymmetry" FP. In addition for U > the line 
of "decoupled chain" FPs is stable. This scenario holds 
up to = 7r/2, at which the "perfect junction" FPs are 
at g = ±\/3, .9 = 0, and the "maximal asymmetry" FP 
is at g — i. For <\> = ir/2 and U > there is a single 
trajectory that runs along the imaginary axis to infinity 
(leading to G — 0) and does not bend back towards the 
real axis as all other trajectories at U > do. In the 
mapping of the complex plane onto the Riemann sphere 
the g = oo FP (north pole) is part of the projected line 
of "decoupled chain" FPs and shows the same stability 
properties and scaling dimension. 

To obtain the scaling dimensions of the FPs we gener- 
ically analyze the scaling of G v y as a function of Sn 
with respect to the FP conductance which is 1, 4/9, or 
depending on the FP studied. For sufficiently large 
N (in some cases up to 10 5 sites are required) we find 
power-law scaling close to all FPs. In cases were the FPs 
are not stable this scaling does not present the asymp- 
totic behavior, but exponents smaller than can still 
be read off at intermediate N. Approaching the "per- 
fect junction" FPs along C(<p) for > Eq. Q yields 
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G v y — 4/9 cx Img, i.e. G v . v < — 4/9 and Img scale with 
the same exponent. The scaling dimensions of the "per- 
fect junction" FPs at <j> = cannot be read off from the 
conductance as G = 4/9 on C(0) and we use the scaling 
of Img. 
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FIG. 3: Scaling exponents close to FPs. 

The scaling exponent of the "maximal asymmetry" FP 
is independent of the flux <f> > and the direction from 
which it is approached (U < 0) or in which it is left 
(U > 0). Its [/-dependence is shown in Fig. |3l (squares) 
and agrees to leadin g or der in U with the prediction 
2(A - 1) (solid Une),|l2] with A = AK/(3 + K 2 ) and 
K given in Eq. (J2J) . In our scheme terms of order U 2 are 
only partly included and we cannot expect agreement to 
higher order. The non-monotonic behavior of 2(A— 1) is 
not reproduced and in our approximation the "maximal 
asymmetry" FP stays attractive for all U < 0. For small 
to intermediate \ U\ with TLL parameters 1/2 < K < 3/2 
we confirm the conjecture of Ref. 

For the line of "decoupled chain" FPs we find for all <j> 
the same scaling exponent j3 s , as found for the single im- 
purity problem close to the respective "decoupled chain" 
FP, applying the fRG. Its dependence on U is shown in 
Ref. |a and agrees to leading order with 2{1/K — 1). 

As discussed above the = "perfect junction" FP 
at g = 1 [yellow diamond in Fig. |21 a)] has properties 
different from those of the other "perfect junction" FPs 
(blue circles in Fig. The scaling exponent 7 of the 
latter FPs is independent of 4> and shown in Fig. [31 (cir- 
cles). To leading order it is given by 7 w U j (37r) (dashed 
line). This form does not coincide with the expansion 
of any of the above X-dependent exponents [after using 
K w 1 — U/n; see Eq. (J2J]. We thus find a new scaling 
dimension. Since higher order terms are only partly in- 
cluded we cannot determine the functional dependence 
of 7 on K . To derive such an expression presents a chal- 
lenge for methods that do not require approximations in 
the strength of the interaction. For the cf> — "perfect 
junction" FP at g — 1 the scaling exponent (of Im g) is, 
up to our numerical accuracy, equal to —7. For junc- 
tion parameters that initially do not fall onto C(0), but 
are close to it, one can read off an exponent from the 
conductance that describes how C(0), with G = 4/9, is 



approached (U < 0) or left (U > 0). It is equal to the 
fRG exponent j3 w of the "perfect chain" FP in the single 
impurity problem. Its [/-dependence is presented in Ref. 
and agrees with 2(K — 1) to leading order. 

The appearance of stable FPs with symmetric con- 
ductance G — 4/9 at U > Oisa surprising result. 
Even though time reversal symmetry is explicitly bro- 
ken at 4> > 0, due to correlations the conductance of 
systems with parameters on C(<j)) behaves as if time re- 
versal symmetry is restored. It is remarkable that close to 
the line of "decoupled chain" FPs the relative difference 
\l{G v y + G„y) scales as Spf and thus 
vanishes if U > 0. This implies that G u y and G„v be- 
come equal faster than they go to zero. In that sense also 
on this line of FPs and thus for all junction parameters 
time reversal symmetry is restored if U > 0. 

Using the fRG we determined the complete renormal- 
ization group flow for a TLL Y-junction pierced by a mag- 
netic flux. Besides uncovering a new type of low-energy 
physics we demonstrated the power of our approximation 
scheme. Usually field theoretical models are used to in- 
vestigate the transport properties of inhomogeneous TLL 
applying methods that are specific to such models. This 
way exact results for either fairly simple geometries (sin- 
gle impurity 0|) or restricted parameter regimesflij were 
obtained. Our method can be applied to microscopic 
models with arbitrary junction parameters and provides 
results for the conductance that are accurate for small to 
intermediate \U\. It can also be used to study transport 
on intermediate and large energy scales. [S| Furthermore, 
the technique can be applied to investigate the transport 
through more complex networks of TLLs that might be- 
come important in future nano-electronic applications. 
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* Xavier B.-T. passed away in a tragic traffic accident on 
August 15, 2004. 
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